Quantum nucleation in a single-chain magnet 
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The field sweep rate (v = dH/dt) and temperature (T) dependence of the magnetization reversal of 
a single-chain magnet (SCM) is studied at low temperatures. As expected for a thermally activated 
process, the nucleation field (H n ) increases with decreasing T and increasing v. The set of H n (T, v) 
data is analyzed with a model of thermally activated nucleation of magnetization reversal. Below 
1 K, H n becomes temperature independent but remains strongly sweep rate dependent. In this 
temperature range, the reversal of the magnetization is induced by a quantum nucleation of a 
domain wall that then propagates due to the applied field. 

PACS numbers: 75.10.Pq, 75.40.Gb, 76.90.+d, 75.45.+j 



Recent efforts in synthetic chemistry has led to a 
quickly growing number of magnetic systems that show 
slow relaxation of magnetization. Apart of interests for 
applications in, for instant, ultrahigh density magnetic 
recording, such systems are ideal to test theories. A well- 
known example is the single-molecule magnet (SMM) 
that exhibit slow magnetization relaxation of their spin 
ground state, which is split by axial zero-field split- 
ting [l|, [2( . SMMs are among the most promising candi- 
dates for observing the limits between classical and quan- 
tum physics. A more recent example is the single-chain 
mag net (SCM) HH0 showing slow relaxation of mag- 
netization as the consequence of the uniaxial anisotropy 
seen by each spin on the chain and magnetic correla- 
tions between spins. Although it seemed that there was 
a reasonable agreement between the experimental data 
and Glauber's theory of a ID Ising spin chain 0], it 
was shown that several other arguments should be con- 
sidered to fill the gap between the theory and the ex- 
perimental results [7J. The most important arguments 
concerned the introduction of magnetic anisotropy and 
finite-size effects. Indeed, their influence on the static 
and dynamic properties of the SCMs was confirmed ex- 
perimentally 0,SH[I3- 

Quantum tunneling of domain walls in a ID meso- 
scopic ferroma gne tic s amp le was theoretically investi- 
gated 0, 0, 13 E-l E3 an d crossover tempera- 
tures between the classical and quantum regime were 
predicted. Domain wall nucleation and depinning were 
studied in single Ni wires and showed indeed a flatten- 
ing of the temperature dependence of the mean switching 
field (H sw ) below about 5 K J7\ and 1 K 0. Because 
of surface roughness and oxidation, the domain walls of 
a single wire were trapped at pinning centers. The pin- 
ning barrier decreases with an increase of the magnetic 
field. When the barrier is sufficiently small, thermally 
activated escape of the wall occurs. This is a stochastic 
process that can be characterized by a switching (depin- 



ning) field distribution. A flattening of the temperature 
dependence of H sw and a saturation of the width of the 
switching field distribution were observed. The authors 
proposed that a domain wall escapes from its pinning site 
by thermal activation at high temperatures and by quan- 
tum tunneling below T c ~ 5 K and 1 K How- 
ever, such crossover temperatures are about three orders 
of magnitude higher than the T c predicted by current the- 
ories |15J. The propagation of a domain wall across an 
energy barrier in a domain wall junction was also stud- 
ied and preliminary investigations seems to indicate the 
possibility of quantum tunneling below 0.7 K 

In this letter, we show that tunneling can occur in 
truly ID systems like SCM provided that a driving field 
is applied that lowers the energy barrier. Indeed in zero 
applied field, the probability of tunneling is negligeable 
due to the exponential increase of the correlation length. 

The studied SCM is a heterometal- 
lic chain of Mn 777 and Ni 77 metal 
ions : [Mn 2 (saltmen) 2 Ni(pao) 2 (P3O2] (CIO4 ) 2 

(saltmen 2- = N,N'-(l,l,2,2-tetramethylethylene) 
bis(salicylideneiminate); pao - = pyridine- 2-aldoximate; 
Py = pyridine), called M112N1 chain henceforth 0, 0. 
At low temperatures, this compound can be described 
as a chain of ferromagnetic coupled 5 = 3 [Mn 777 -Ni 77 - 



Mn 



mi 



units. 



The spin system of the M^Ni chain can be described 
by an anisotropic Heisenberg model: 

H = - J 22 SiS i+1 - D \] S? z - guBlio^SiH (1) 

i i i 

where J is the ferromagnetic exchange constant between 
the spin units and D is the single-ion anisotropy. D/fce 
= 2.5 K was obtained from magnetization measurements 
as a function of a magnetic field applied perpendicular to 
the easy axis 0. AC and DC relaxation time measure- 
ments showed a unique relaxation time over 10 decades. 
Above 2.7 K, the thermal dependence of the relaxation 
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time followed an Arrhenius law with an activation energy 
of 74 K. Below 2.7 K, a departure from this simple be- 
havior was observed and a smaller activation energy of 
55 K was found around 2 K 0. The crossover at about 
2.7 K was interpreted as the manifestation of finite-size 
effects 0- Indeed, the activation energy of the relaxation 
time should decrease from (4J + D)S 2 to (2 J + D)S 2 at 
the temperature where the correlation length equals the 
chain length 7]. The exchange energy J/k-Q = 1.56K was 
found which was in agreement with independent thermo- 
dynamical measurements. A saturation of a semi-log plot 
of \T versus 1 /T was used to estimate the mean chain 
length of about 100 Mi^Ni units. 

The relaxation rate at H = is extremely small be- 
low 1.4 K. We applied therefore a magnetic field to study 
the low temperature relaxation process. The magnetiza- 
tion measurements were performed by using (i) a magne- 
tometer consisting of several 10 x lOum 2 Hall-bars |20| 
and (ii) an array of micro- SQUIDs |2J on top of which 
a single crystal of Mr^Ni was placed, for higher and 
lower fields than 1.4 T, respectively. The field can be 
applied in any direction by separately driving three or- 
thogonal superconducting coils. The field was aligned 
with the easy axis of magnetization using the transverse 
field method [12 . 

Typical hysteresis loops are presented in Fig. ^ The 
Mn2Ni chain displays smooth hysteresis loops which are 
strongly temperature and field sweep rate dependent. 
The temperature and field sweep rate dependences of the 
coercive fields (called mean nucleation fields H n hence- 
forth) were measured and plotted in Fig El As expected 
for a thermally activated process, H n increases with de- 
creasing temperature T and increasing field sweep rate 
v = dH z /dt. Furthermore, all our measurements showed 
an almost logarithmic dependence of H n on the field 
sweep rate (Fig. Et). H n becomes temperature indepen- 
dent below about 0.5 K. 

We analyzed the set of H n (T,v) data with a model 
of thermally activated nucleation of magnetization rever- 
sal analogous to that of a magnetic single-domain parti- 
cle [H [24], IH |2(| . AC and DC relaxation measurements 
at H = showed that the magnetization reversal below 
2.7 K is dominated by the ends of the Mr^Ni chain Q and 
the dynamics were described by the generalized Glauber 
model taking into account chains with finite lengths. At 
sufficiently low temperatures and at zero field, the energy 
barrier between the two states of opposite magnetization 
is much too high to observe a reversal. However, the bar- 
rier can be lowered by applying a magnetic field in the 
opposite direction to that of the chain's magnetization. 
When the applied field is close enough to the nucleation 
field of a domain wall, thermal fluctuations are sufficient 
to allow the system to overcome the nucleation barrier, 
and a domain wall nucleates. Then, due to the applied 
field, the magnetization of the entire chain reverse via 
a domain wall propagation process. The domain wall 
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FIG. 1: (color online) Hysteresis loops for the Mn2Ni chain at 
(a) several temperatures and 0.008 T/s and (b) at 0.04 and 
1.4 K and at several field sweep rates. 



nucleation can be thermally activated at high tempera- 
tures or driven by quan tum tunneling at low tempera- 
tures m ei in mm m 

This stochastic nucleation process can be studied via 
the relaxation time method consisting of the measure- 
ment of the probability that the magnetization has not 
reversed after a certain time. In the case of an assembly 
of identical and isolated spin chains, it corresponds to 
measurements of the time dependence of magnetization. 
The probability that the magnetization has not reversed 
after a time t is given by: 

P(t) = e^' T (2) 

and t can be expressed by an Arrhenius law of the form: 

t{T,H) = To e AE ^^ T (3) 

where AE(H) is the field dependent nucleation energy 
barrier and tq is a prefactor which is supposed to be a 
constant. In most cases AE(H) can be approximated by: 

AE(H) a E (1 - h) a (4) 
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FIG. 2: (color online) The mean nucleation field H n for the 
Mn2Ni chain as a function of (a) temperature and (b) field 
sweep rate. 



where h = H/H®, H® is the nucleation field at zero tem- 
perature, Eq is roughly the nucleation barrier height at 
zero applied field, and a is a constant of the order of 
unity (for most cases 1.5 < a < 2 

There are two limiting cases: (i) for D << J, the 
nucleation energy barrier can be calculated exactly |2rl 
MM: 



AE(H) = 2/cV2J£>(tanLR - hR) 



(5) 



where R — arcosh(^/l/ft.) and k — 1 or 2 for the nu- 
cleation at the ends or insight the chain, respectively. 
For h = Q, this energy is the well-known energy of 
one or two domain walls, respectively. For h — > 1, 
AE(H) w 4fc/3 v / 2JD(l - hf/ 2 - (ii) for J = 0, the spins 
are decoupled and we have the case of a Stoner-Wohlfarth 
particle |3JJ, |31| with uniaxial anisotropy. When the field 
is applied along the easy axis of magnetization, all con- 
stants can be determined analytically jU a = 2, 
E Q = KV, and the switching field H° w = 2K/M S , where 
K is the uniaxial anisotropy constant, V is the par- 



ticle volume, and M s is the saturation magnetization. 
For SMMs with dominating uniaxial anisotropy: a = 2, 
E Q = DS 2 , and iJ s ° w = 2DS/ gt i of , B . 

In our case of the Mn2Ni chain with D/k B = 2.5 K and 
J/k B = 1.56 K we are not aware of an analytical expres- 
sion and propose to use a = 2 because of the arguments 
developped in [29| . 

In order to study the field dependence of the relax- 
ation time t(T, H) and to obtain the parameters of the 
model, the decay of magnetization has to be studied at 
many applied fields H and temperatures T. This is ex- 
perimentally very time consuming. A more convenient 
method for studying the magnetization decay is by ramp- 
ing the applied field at a given rate and measuring the 
mean nucleation field H n which is the field value to obtain 
zero magnetization (coercive field) . H n is then measured 
as a function of the field sweep rate and temperature 
(Fig |2J . An analogous procedure 0,|U|3j| was applied 
to nanoparticles [34| and recently to SMMs [35|. The 
mean nucleation field of an assembly of identical non- 
interacting SCMs is given by: 



H n (T,v)^H°[l- 



En \v 



l/a N 



(0) 



where the field sweeping rate is given by v = 
dH z /dt; H® is the nucleation field at zero tempera- 
ture, and c depends on the details of the approxima- 
tions: c = H°k B T/[T aE (l - H.JH^- 1 ] in 0, 
d = HKEo/kTf^/iTQa) in [13, and it can be taken 
constant when the exact value of To is not needed. We 
applied the three approximations to nanoparticles [34| 
and to SMMs |3^| and found that the first approxima- 
tion gives a To which is closest to that extracted from an 
Arrhcnius plot. 

The validity of Eq. [S] was tested by plotting the set 
of H n (T,v) values as a function of [Tln^/v)] 1 ^ 2 where 
c = H°k B T/T a 2E a (l - H n /H°). If the underlying model 
is sufficient, all points should collapse onto one straight 
line by choosing the proper values for the constant To- 
We found that the data of H n (T, v) with T > IK fell on 
a master curve provided To = 7.4 x 10 -9 s (Fig.[3|). 

At lower temperatures, strong deviation from the mas- 
ter curves are observed. In order to investigate the pos- 
sibility that these low-temperature deviations are due to 
escape from the metastable potential well by tunneling, 
a common method for classical models is to replace the 
real temperature T by an effective temperature T* (T) in 
order to restore the scaling plot 0, |35| . In the case of 
tunneling, T*(T) should saturate at low temperatures. 
Indeed, the ansatz of T*(T), as shown in the inset of 
Fig. 03, can restore unequivocally the scaling plot demon- 
strated by a straight master curve (Fig. Oh)- The flat- 
tening of T* corresponds to a saturation of the escape 
rate, which is a necessary signature of tunneling. The 
crossover temperature T c can be defined as the temper- 
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[T*ln(c/v)]0.5 m,H n (T) 



FIG. 3: (color online) (a) Scaling plot of the mean nucleation 
field H n (T, v) for field sweep rates between 0.001 and 0.14 T/s 
and several temperatures: 0.04 K and from 0.2 to 2.6 K in 
steps of 0.1 K. (b) Same data of H n (T, v) and same scales but 
the real temperature T is replaced by an effective temperature 
T* (see inset) which restores the scaling below 1.1 K. 



FIG. 4: (color online) (a) Field dependence of the energy 
barrier of the Mn2Ni chain obtained from Eg. 1101 and the set 
of Hc(T, v) data from Fig. El (b) Same data of H C (T, v) but 
the real temperature T is replaced by an effective temperature 
T* (see inset of Fig.Eb). 



ature where the quantum rate equals the thermal one. 
The inset of Fig. Ob gives T c = 0.8 K. The slope and the 
intercept of the master curves give E^/k-Q — 47 K and 
fM)H° = 6.95 T. 

Several points should be mentioned: (i) Eq. is not 
valid for fields which are close to H = because the 
model only takes into account the nucleation from the 
metastablc to the stable state. However, close to H = 0, 
transitions between both states are possible leading to a 
rounding of the master curve at small fields; (ii) the re- 
maining energy barrier for the tunnel process at 0.04 K 
can be calculated using Eq.0] We found ~30 K for a field 
sweep rate of 0.001 T/s; (iii) in case of a distribution of 
nucleation barriers, different parts of the distribution can 
be probed by applying the method at different M values; 
(iv) this method is insensitive to small intermolecular 
interactions when H n is much larger than the typical in- 
teraction field; and (v) the method can be generalized for 
2D and 3D networks of spins. 

In conclusion, the presented low-temperature studies 



of the field driven magnetization reversal of the Mn2Ni 
SCM suggest for the low temperature region that the 
magnetization reversal starts by a quantum nucleation 
of a domain wall followed by domain wall propagation 
and reversal of the magnetization. Further studies will 
concern the application of transverse fields which should 
enhance the quantum nucleation rate. 

This work was supported by the EC-TMR Network 
QuEMolNa (MRTN-CT-2003-504880), CNRS, Rhone- 
Alpes funding, the Universite Bordeaux 1, and the Con- 
seil Regional dAquitaine. 



ANNEX 

The probability density of reversal of a stochastic pro- 
cess is 
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and the maximum of the probability density can be de- 
rived from 

d 2 P P ( dr\ 

This gives the following general result for the maximum 
of the probability density 



The application of the result to Eq. leads to 

AE(H) = fc B Tln ( -^F j (10) 

Using Eas.HlandlTni we find Eq. El 

Eq.E3can be used to plot directly the field dependence 
of the energy barrier (Fig. 0} 
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